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ABoITRACT 


A technique of adjusting the initial steering currents to improve 
the short-term forecasts of tropical cyclone motion with a dynamical 
model was tested. Initial tests were performed with hand-analyzed data, 
although 41 objectively-analyzed cases were also examined to evaluate 
Operational feasibility. Several schemes to track the storm objectively 
resulted in the use of the streamfunction minimum to define the storm 
centers. 

Two blasing approaches were developed. The first used statistical 
regression equations to predict the initial storm motion, while the 
second used an empirical relationship to define the initial steering 
current. Various areal distributions of the bias corrections were then 
applied to the initial wind fields. Results from these tests indicated 
improved forecast accuracy if the largest corrections were applied at 
some distance fram the storm center. The biasing decreased average 
forecast errors for the first 30 hours by about 40% for the hand-analyzed 
cases and about 30% for the objectively-analyzed cases, with 25% and 
15% decreases after 48 hours, respectively. Significant improvements 
were also observed in the standard deviation of track errors, with 
averaged decreases of 60% for hand-analyzed cases and 45% for objectively- 


analyzed cases. 





iL lee 


iid, 


Jig A 


ti-sbp Or CONTENTS 


INTRODUCTION - - - - ----fee87F724-842 22 2S 
HE TROPICAL CYCLONE MODEL = =---+-+-+-+<-s- - - -18 
A. THE MODEL - =----+-+efffe7f+277+2e 18 
B. MODEL MODIFICATIONS - - ----+e-+-+e-e-e- - 18 
C. THE TRACKING METHOD - - - - =-+ee-e-e-e-- 20 
INITIALIZATION - - -- ---+e-ee-e-eeFFf-+2e2-° Ze 
A. TROPICAL VS EXTRATROPICAL - - - ---+-e-- oy 
B. DATA INPUT - - - - - ------- - ee 
Daeveher ten) OF THE BIASING PROCEDURE - - =- = = = 23 
A. THE PROBLEM - - - - = =- = --+-+e-ee-e--e- - fase 
B. THE BIAS INPUT - - - - -----+----e-- aS, 
C. DEFINITION OF THE STEERING CURRENT - - - - - 30 
Pee oe Me PeONS OE THE EXPECTED 
MODEL STORM MOTION - - - - = ----+-+-e- - Sh, 
1. Initial Observations - - -------- 30 
2. Model Characteristics -----+---- ait 
3. The Prediction Equations - - - --- - - 41 
E. THE BIAS CORRECTORS - - - - =--+-+-+--+e-- 43 
1. The First Approach - - ---+---+-+--- 4 8 
2. Second Approach - ---+-+e-+-+-+-+-- - 48 
F. THE BIASING SCHEME - - - - --+--+-e-+-e-+e--- - 20 
1. Initial Experiments - -----+-+--- 50 
2. Differential Adjustments - - - - -- - - 29 





V. RESULTS - ----7+7727777 77 - - - ee - 57 


RIC AGESG S=)- = == - = = = = eee 
Phere AGECm=e = 22 = = = & = = - = - -~ 63 
‘tT, 8 «OCWSTnG = 3.5 30S 78 
eee pe THe TROPICAL CYCLONE MODEL - - - - -- - = - 81 
&, ME MODELM G¢a@e Se 81 
2, SS. 62D Sele ee 82 
Pee OU RCONDITIONS =.= = =.- - =< SAEs 2 83 
Pe ENC SmmROencAl CYCLONE MODEL’ - - -- = - - 85 
PIC OSMEONEHASH = a0] 2 ~ = = - Sv- = Se = 86 
DAMPER EEERENORS 2 = - = = = === === --= = + == 90 
DUPE DESOneTBUTTON LIST - - —------+-+--- 92 





JS 


Pio Ob ABLES 


Mean 6-hourly forecast errors (km) (I) based 

on Best track for the 1974 cases for the 

unbiased, first bias and second bias approaches. 

The error at each 6-h interval is normalized by 

the initial position error. Also tabulated 

are the number of cases (II) which equal or 

exceed the unbiased error - - -----+-+-+-+- - 64 


same as Table I except for the 1975 cases and 
without the second bias method - - ---+--+-+---- 73 





IO 


a 


Pitot or FIcuRES 


Comparison of forecast tracks with and 
without vortex heating starting from 
00 GMT on 25 Nov. 1974 - - = = - = - - 


Predicted 850-mb wind field at 30 h 
starting from 00 GMT 27 Nov. 1974 with 
vortex center determination from vor- 
ticity, geopotential, streamfunction 
and enhanced streamfunction fields = - 


Storm tracks based on vorticity, geo- 
Perenttadileancd StreamtUnetlon fileids for 
iS oncem (ric. 2/ Nov, I19/t - - — = = = 


Adjusted wind components at near-zero 
wind speeds derived from the biased 
inverse enhancement function (see text) 


Average u and v components over succes- 
Sively larger domains for three levels 
based on the 1974 cases - = - - - = = 


same as Fig. 5 except for '1975" - - = 


Average initial steering currents, best 
~mack and model storm motion in first 
Six hours for the 1974 cases - - - = - 


Same as Fig. 7 except for '1975' - - - 


Empirical mean flow relationship with 
Bespeen tO actual Storm motion vector 
based on George and Gray (1976). Mean 
flow is 86% of 700-mb speed and 16” to 
the right of 500-mb direction - - - - 


850-mb wind field for Typhoon Gilda, 
O04 August 1974 ----+-+-+-+--+-- - 


Non-divergent winds at 850-mb for same 
Field as in Fig. 10 ---+-+-+-+-4-- 


Same as Fig. 7 except with steering 
currents based on non-divergent winds 


ae 


oak 


Ze 


Zo 


32 


34 


So 


Sy 


4 Q 


42 


43 


uy 





se 


ee 


Ise 


ie. 


ey. 


ne). 


iS. 


ZO). 


alle. 


Lage 


Zee 


24. 


LS. 


Zs 


Same as Fig. 8 except with non-divergent winds - 


Senematwec Or ene bias corrector definition 


from (a) the first approach using the predicted 


initial storm motion and (b) second approach 
using the empirical steering current - - - - 


Examples of various horizontal distributions 
of the bias corrector with respect to the 
grid point nearest the storm center. In 
each case a unit value corresponds to the 
desired bias corrector as in Fig. 14 - - - - 


Same as Fig. 1l after application of the 
betas “COrrPeectoo = = = = = = =- = = = = = a ages 


Model forecast tracks from 00 GMT O04 Aug 
1974 versus (A) best track and (B) unbiased 
ema fOr em (CG) ) box-like, (Dp) Cosine and 

2) baste Dids fUnetion as in Fig. 15e€ - - = 


Biasing function for a v-component bias 
corrector derived after differential adjust- 
ment tests with the 1974 cases - ----- - 


Model Saeorecasts with bias distribution as 
in Fig. l5e for the 01-04 Aug. 1974 cases - 


Average forecast error (km) with 1974 cases 
for JTWC (Ce), unbiased model (a), first 
bias (#*) and second bias (* ) techniques. 
The error at each 6-h interval is normalized 
by the initial position error - - - --- - 


Same as Fig. 20 except for standard devia- 
tion of track errors (km) ---+-+---- = 


Same as Fig. 18 except for 1975 cases, 
including the positive and negative zonal 
S®MpOneneTeCOrrectcors = = =—=- == - =. --- 
Peak biasing values for the 1975 cases as 

a function of direction and magnitude of 
Pepto meckousm (isa m- = = = - = = == 
same as Fig. 20 except for 1975 cases - - - 


same as Fig. 21 except for 1975 cases - - - 


same as Fig. 24 without normalization by 
initial position error - - - - - - Bee eh pe 


Vertical distribution of dependent variables 


and pressure levels for the three-dimensional 
model (after Harrison, 1973) -----+-+----- 


US 


Ug 


o1 


54 


98 


60 


62 


65 


66 


69 


7k 
74 


(is 


76 





eo 


| 


Pilot Oreo IbOLS 


Time variable 

Velocity component in the x-direction 
VEUGetEy geOmponen. In the y-direction 
West-e€ase, MNOrth=sSsouth space axes 
Potential temperature 

Coriolis parameter 


Components of shear stress per unit area by 
u Wind component 


Components of shear stress per unit area by 
Vv Wind component 


Acceleration of gravity 

Map factor 

Vertical distance above 1000 millibar surface 
Geopotential height 

Geopotential height of 1000 millibar surface 
Vertical velocity in pressure coordinates 
Pressure (space axis in vertical) 
Seneanmrhuncrlen 

VOeCeriLca Gy 

Specific heat at constant pressure 

Finite difference operator 

Gas constant for dry air 

Atmospheric density 


Average of S over the entire grid 


LO 





Av 


ei 


<| 


Average of S along one row of the grid 
Perturbation component in the x-direction 
Laplacian operator 

Predicted model component in the x-direction 
Predicted model component in the y-direction 
Dw otcasceceorrectOr in the x-direction 

Unit bias corrector in the y-direction 


Vertically averaged steering mean current 
component in the x-direction 


Vertically averaged steering mean current 
component in the y-direction 


ak 





ACKNOWLEDGEMENTS 


The author wishes especially to thank Professor R. L. 
Elsberry for his patience and guidance throughout the course 
of this research. 

Special appreciation 1s extended to Lt. D. E. Hinsman, 
USN, whose concerned interest and continual assistance en- 
abled the majority of the tests and results presented herein 
Bessible. 

The author also wishes to thank Professors G. J. Haltiner 
aag@ec. FF. Chang for reviewing the manuscript and making 
several useful suggestions. 

The personalized service received from the W. R. Church 
Computer Center night shift personnel gave invaluable numeri- 
Smeesupport. Lts. (USN) W. F. Johnson and O. M. Lubeck 
extended expert computer assistance in the solution of 
Several technical difficulties. 

An affectionate thank you is extended to the author's 
wife, Judy, for her devoted patience during the research and 
unyielding assistance with preliminary typing of the 


Manuscript. 


eZ 





Pe LN ORC Teoh 


Prediction of tropical storms 1S a major concern world- 
wide. The need to know storm conditions in advance has pro- 
duced a variety of prediction methods. Subjective methods, 
including persistence, have normally produced good short- 
range forecasts. The occurrence of storm recurvature 1S an 
especially difficult problem for any prognostic scheme, 
whether subjective or objective. This was a significant 
problem encountered during the 1975 western Pacific typhoon 
season. An unusually large number of the storms tracked 
north, while some quickly recurved. As a result, all error 
statistics were higher than normal. Meanwhile, much progress 
has been made in the numerical simulation of characteristics 
of tropical storms, such aS intensity and even spiral bands 
(Anthes et al, 1974), but the major emphasis has focused on 
forecasting storm movement (Hovermale et al, 1976, and Ley, 
1975). A prime objective of the latest multi-level, nested- 
grid models has been to improve medium- and long-range motion 
forecasts over those of persistence, analog, or statistical 
forecast techniques. Yet, no matter how sophisticated the 
numerical models may be, initialization with limited and poor 
data may produce undesirable consequences. Even if the 
analysis is good, the model should at least be able to simu- 


late the initial storm motion with reasonable accuracy. It 


ire 





is observed that dynamical models occasionally produce ir- 
regular storm tracks, and the forecast motion is generally 
too slow. A model cannot be expected to operate at its full 
potential and produce extended, accurate forecasts if the 
short-term forecast is bad. Specification of the initial 
data fields seems to be a likely source of error in these 
Seeuations. 

Since the resources required to increase the amount of 
high quality data are costly and will require considerable 
amounts of time to acquire, alternative methods for improv- 
ie. initialization of models should be considered. Of the 
resources currently available, the valuable information re- 
presented by the recent movement of the tropical storm has 
Wem to be effectively utilized in a dynamical model. All 
meristical schemes for predicting tropical storm motion 
(Renard et al, 1972) rely heavily on the current movement of 
the storm. Perhaps the most direct approach for improving a 
Zmepieal cyclone model is to correct, or bias, the initial 
wind fields to reflect the recently observed motion of the 
storm. It 1s therefore the primary objective of this thesis 
to develop and evaluate a method of biasing the initial data 
to improve the tropical cyclone model forecasts. 

The primitive-equation, three-layer, tropical cyclone 
model (TCM) developed by Elsberry and Harrison (1971) and 
Harrison (1973) was used with a uniform coarse grid. Although 
this model is capable of triply-nested operation (Harrison, 


1973), results from Ley and Elsberry (1976) show that the 
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coarse and nested grids produced nearly identical results in 
a selected case study based on hand-analyzed data. Because 

of the new biasing scheme tested, the simpler, uniform grid 

model allowed for easier modifications and less computer re- 
quirements. 

A bias corrector is desired that would adjust the initial 
data fields initially to steer the storm according to its 
recent movement. If one could estimate the initial storm 
displacement as part of the diagnostic phase of the model, 

a comparison could then be made with the actual storm move- 
ment. If, upon comparison, the storm's predicted initial 
motion deviated significantly from the recently observed 
motion, a bias could be applied to correct the initial model- 
predicted storm track and the subsequent forecasts. In this 
research, the initial model storm movement is calculated from 
a statistical regression equation, and the bias is applied to 
the initial wind fields. To accomplish this, a steering 
current relating the storm movement to the surrounding wind 
Fields was first defined. Then, based on prior storm track 
information, a regression equation was developed to predict 
the storm movement from the steering currents in the three- 
layer model. After comparison of the observed and predicted 
storm motion and determination of the required bias, some 
experimentation was necessary to derive the form of the bias 
tO De applied to the original data fields. In this way, the 
adjusted steering current could steer the storm in the 


desired direction. 
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The biasing method developed was tested on 48 cases in 
the western Pacific region, most of which were of typhoon 
intensity. Seven cases (two storms) used hand-analyzed data 
supplied by the Joint Typhoon Warning Center (JTWC), Guam, 
while the other 41 cases (13 storms) were based on opera- 
tionally-analyzed data obtained from the Fleet Numerical 
Weather Central (FNWC), Monterey, California. Since the 
hand-analzyed cases provided more accurate initial wind 
fields as compared to FNWC's objective analyses (Elsberry, 
1977), they were used for the majority of the tests of the 
biasing method. The contrasting tracks of these two 1974 
storms (one moved westward while the other recurved) also 
provided the directional diversity needed to develop a 
flexible biasing scheme. But the most important result re- 
quired of this scheme was the need to improve consistently 
time EPCM's initial forecasts. Upon completion of the testing 
of the various biasing techniques using the hand-analyzed 
cases, further tests with the operationally-analyzed cases 
from 1975 were evaluated for operational feasibility. All 
forecast positions, biased and unbiased, were verified 
against the best tracks from the Annual Typhoon Reports 
(1974, 1975) published by JTWC. 

Two approaches for biasing the initial data to improve 
the steering current were used. The first approach, which 
was previously described, compared an observed with a pre- 
dicted storm track to determine the required bias. The 


second approach incorporated an empirical relation between 
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the storm motion and the steering current based on the re- 
Sults by George and Gray (1976). As a result of these rela- 
tionships, the empirical steering current was known, and 
could be compared to the observed current to produce a bias 
corrector. The methods and results from these two approaches 


are presented for comparison. 
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ieee Kor ben cy CLONE MODEL 


eee LHe MODEL 

The Tropical Cyclone Model (TCM) characteristics are 
presented in the Appendix. A more complete discussion may 
be found ae Ley (1975). The TCM is a three-layer, primitive 
equation model in pressure coordinates. It is a channeled 
model with free slip conditions on the north and south walls 


and cyclic conditions on the east and west boundaries. 


Eee MODEL MODIFICATIONS 

Since this model was primarily developed for movement 
forecasts, a crude, five-grid point, heating function was 
Semeered over the storm to counteract the dispersion of the 
vortex due to the finite differencing (Ley and Elsberry, 
1976). Because the heating pattern was centered on the grid 
point nearest the pressure minimum, it was sometimes applied 
to the rear quadrants of the storm's direction. Consequent- 
ly, the track of the center tended to zigzag or oscillate. 
A comparison of the 48-hour tracks of one of the storms with 
aaeewiltnout heating (both without friction) 1s given in 
Fig. 1. In almost all cases, this heating caused a retarda- 
tion of storm speed and directional backing with respect to 
the non-heated cases. Therefore, all subsequent case studies 


were computed without the effects of heating or friction. 
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WITH HEATING »s——a 
WITHOUT HEATING e——e« 


NOUNS One 


Fig. 1. Comparison of forecast tracks with and 
without vortex heating starting from 
VOmMeMewr em 2s NOV loys, 
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C. THE TRACKING METHOD 

An objective method of tracking the storm's movement was 
needed for consistent comparisons. A tracking algorithm, 
developed by Lt. D. Hinsman, USN, at FNWC, interpolates the 
maximum and minimum values of any data field. Sequences of 
Six-hourly, model-predicted vorticity, geopotential and 
streamfunction fields defined by 

Seopa VOrGtedey (6) = av/ax = dU/ dy Ci) 

1000 mb, Geopotential (6) = gz (2) 

850 mb, Streamfunction (yp), 

where vy Gs] Ov) Ox — moll 7 oy Ee, 
were used to estimate the storm center movement (Fig. 2). 
Fig. 3 compares the 48-hour center forecasts by the three 
methods tested. The initial storm center was located 
accurately in all fields, since the storm was initially 
bogused symmetrically (reference Numerical Environment 
Eme@edauets Manual, 1975). 

Each method of storm-center tracking was subjectively 
evaluated. After an approximately 18-hour forecast period, 
Significant asymmetrical troughing sometimes appeared in the 
wind fields. Since the coarse-mesh grid was unable to 
resolve the storm's radius of maximum winds, these inner 
winds may become quite weak. The shear and curvature contri- 
butions to the vorticity were much greater in the induced 
Wind troughs. As a result, the vorticity centers may become 
greater in the trough as far as 600 km from the storm's wind 
circulation center. Hovermale et al (1976) experienced this 


tracking problem even in a model with 60-km grid spacing. 
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TYPHOON IRMA - 27 NOVEMBER 1974 

INITIAL POSITION x 

6-HR FORECAST POSITIONS: 
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Fig. 3. Storm tracks based on vorticity, geo- 
potential and streamfunction fields for 
Tyenconva rma. 2/ Tov. 19/74. 
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The second method of tracking using the 1000-mb height 
Minimum increased tracking accuracy only slightly. The u and 
v components defining the anomalous vorticity centers were 
reflected in the @ fields through the equations of motion. 
Therefore, secondary pressure minima were also produced in 
mee troughs. 

lime third method of fem cI AS the storm with its associa- 
ted yw center achieved the most consistent estimates of the 
desired wind field center. These were generally within a 
one-half gridlength (approximately 100 km) error limit. Be- 
fore calculating Vw, the forcing function (f) was partially 
Beeeanetion of uw’, the u-component perturbation field. By 
eliminating the mean zonal flow, circulation centers became 
more easily defined. The streamfunction minimum defines a 
wind circulation center, whereas a vorticity center does not 
necessarily define a low center. Of the three methods tested, 
the w fields provided the most accurate tracking. The Sanders 
integrated barotropic (SANBAR) model, currently in use by 
the National Hurricane Center (NHC), Miami, Florida, also 
uses the streamfunction for tracking (Pike, 1972). SANBAR 
identifies a storm center with a local minimum in w or maximum 
in Vu, but in practice, the average position between the two 
extremes 1S used. 

In some cases with wind speeds less than l ms 71, the 
Streamfunction's minimum value was also weakly defined. One 
might then consider a method of strengthening the gradients 


of w by artificially increasing the wind speeds across the 
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center location. A "biased inverse enhancement" method was 
tested which further increased the accuracy of locating the 
Circulation center. This method is based on a negative, 
Matural log function. The greatest enhancement occurs to 


z 


Perids less than 5 ms (Fig. 4). Biasing the winds by a 


faetor of 107? 


(this bias proved effective, although larger 
or smaller factors may be used depending in the degree a 
enhancement desired) further enhances these fields and the 
resulting y fields, with the maximum emphasis on highlighting 
weak highs/lows. However, these regions of very weak winds 
near the storm center are highly transient features within 
the general trough. The accurate tracking provided by this 
method and the limited resolution of the coarse grid TCM may 
result in various oscillations that are not representative 
of the envelope of the storm. Therefore, all results pre- 
sented in this study were tracked from streamfunction fields 
calculated as in (3). Note that the "biased inverse enhance- 
ment" method approaches the limits of interpolating fields 
of finite data. The location of extraneous perturbations may 
not be representative of the true circulation. Since the in- 
creased resolution provided by nested-grid models can 
realistically depict small-scale wind fields, those circula- 
tion centers tracked by the "enhanced" method would be more 
meaningful. 

It 1s of interest to note that because of these results, 
tracking tropical storms with p minima has already been incor- 


pesaced amtO FNWC's TCM in their operational support of JTWC. 
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Fig. 4. Adjusted wind components at near-zero 
wind speeds derived from the biased 
inverse enhancement function (see text), 
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This was made possible through the efforts of Lt. D. Hinsman, 


USN, at FNWC, who evaluated and implemented the method. 
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Pye PROPICAL VS. EXTRATROPICAL 

In general, initialization of tropical models differs 
from their extratropical counterparts. In a mid-latitude 
model, the none is normally derived from a solution of the 
balance equation with an objectively analyzed geopotential 
field (Haltiner, 1971). Such a procedure is not suitable for 
a tropical model since the pressure gradient in the lower 
latitudes is very weak and the pressure field is not well 
defined. Moreover, it can be shown that the mass field pri- 
Mariiy tends to adjust to the wind field in the tropics 
during the early stages of the integration after initiali- 
zation, especially in mesoscale models as found by Anthes 
(1974). Consequently, in the tropics, primary consideration 
must be given to the wind field when performing the initial 
diagnostic phase. 

A complete discussion of the theory and finite differ- 
ence equations used for the TCM initialization can be found 


in the Appendix. 


Eee DATA INPUT 

Wind data for the 1974 cases was supplied by the U. S. 
Fleet Weather Central/Joint Typhoon Warning Center, Guam, 
Mariana Islands (FWC/JTWC). The 850-, 500-, and 250-mb 


Streamlines and isotachs were hand-analyzed for Typhoons 
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Irma and Gilda. These analyses, plus careful extraction of 
the data to the grid, provided a data base superior to the 
objective analyses used by the 1975 semi-operational cases 
(Elsberry, 1977). This combination of accurate hand-analyses 
and the directional variations of Irma (westward track) and 
Gilda (recurving track) provided diversified test samples. 
Not only did these two storms follow widely varied tracks, 
but the TCM forecast both good and poor tracks even though 
their initial data was hand-analyzed. It was therefore 
imperative that the biasing scheme presented should not 
markedly alter a good forecast. 

For convenience, the hand-analyzed, synoptic data at 
500- and 200-mb (800-, 700-, 400-, and 250-mb for 1975 cases) 
were moved to the nearest model predicted level (reference 
Appendix, Fig. A-1). Subsequent reference to these levels 
will be according to the synoptic data fields, not the model's 
predictive levels. 

Data initialization for the 1975 cases was from the FNWC 
global band upper air analyses (GBUA). The bogusing of upper- 
level cyclonic outflow has not yet been incorporated into the 
PMetailiZation procedure of the FNWC TCM. For this reason, 
plus the incomplete assimilation of Defense Meteorological 
Satellite Program (DMSP) data (see Elsberry, 1977), these 
cases had significantly inferior upper-level circulations 


compared to the 1974 cases. 
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heer ewer Or nt Ve tAoING PROCEDURE 


pee LHE PROBLEM 

The problem confronted by this thesis was the consistent 
initial forecast error of FNWC's TCM. In some cases, the 
error was small, but in other cases, the model forecast 
track deviated significantly from the most recent movement 
of the storm. The hypothesis is that the observational data 
base was inadequate in these cases, and adjustments should 
be made to the initial data. 

The purpose of this thesis was to bias objectively the 
TCM's initial wind fields from the storm's recent 12-hour 
movement and to test this procedure on a semi-operational, 
real-time basis. To be discussed are the bias input, which 
represents the recent 12-hour movement centered on the 
observation time; the steering current, which interrelates 
the cyclone/environmental flow fields; the statistical 
regression equation, which predicts the TCM's storm movement; 
and the bias corrector and how it was applied to the TCM's 


Tame fields. 


Emeline BIAS INPUT 

The vector which represents the best estimate of the 
recent storm motion was defined as the bias input. The bias 
input would, on a real-time basis, be a past 12-hour vector 


provided by JTWC or any operational site. The 12-hour vector 
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was chosen to reduce errors produced by sporadic and short- 
term oscillations introduced by various positioning techniques 
P@omesynoptic, aircraft, and satellite data. Since FNWC's 

TCM begins at observation time plus 9 hours, the 12-hour 
vector could be centered about the observation time. For 
example, the bias input for a 0000 GMT observation time would 
begin at the 1800 GMT position from the previous day and end 
at the 0600 GMT position. While developing this biasing 
technique, best-track positions were used to insure a good 
bias input. In an operational mode, this vector would be 
derived from tropical cyclone warning positions, which would 
introduce some additional error. For better accuracy, a 

blas input consisting of a best-track position at observation 
time minus 6 hours and the warning position at observation 
time plus 6 hours could easily be developed and utilized at 
an operational site. If such a continually-updated tracking 
program were developed, the accuracy of the bias input would 
improve over those using only warning positions. It should 

be noted that the input vector is linear and therefore becomes 
less reliable as the true storm velocity accelerates, or 


changes direction rapidly. 


See eer INITION OF THE STEERING CURRENT 

Before any modification or correction could be applied 
to the model, a standard, cyclone/environmental wind field 
relationship had to be established. The proper domain of 
horizontally-averaged wind fields were needed to define a 


representative steering current. This averaged current was 
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needed later to estimate the difference between the actual 
storm track and that predicted by the TCM. 

Before determining the steering current, each storm center 
was first assigned an origin at its nearest grid point. 
Second, area averages over increasingly larger domain sizes 
were calculated to find a consistent and representative wind 
field relating cyclone motion to environmental flow. These 
domain sizes ranged from 3 x 3 to 13 x 13 grid squares 
centered on the storm. To simplify the problem, wind averages 
from each level (850-, 500-, and 250-mb) were added to define 
the vertically-averaged steering current. 

Averaged u and v components as a function of increasing 
domain are depicted in Fig. 5 for the 1974 (hand-analyzed) 
cases. Note that the averages show a weak, positive v compo- 
nent and a strong, negative u component reflecting the general 
WNW movement of tropical cyclones. Averaging over larger 
areas should tend to decrease the effect of the dominant 
memoon Circulation, with the result being the steering 
current. The ideal average would produce a slope rapidly 
decreasing with increasing area indicating a region beyond 
some zone where the component averages become invariant. 

This zone would be equally influenced by the cyclone/environ- 
mental flow fields. This zone is evident from the 850-mb's 

u and v components (Fig. 5) near the 7 x 7 (12 degree) square 
which is approximately equal to a circle with a 6%° radius. 
Jones (1976) also observed this zone at a radius of about 


600 km using a triply-nested model with a middle grid 
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resolution of 30 km. Similar results were obtained by George 
and Gray (1976) by averaging rawinsonde data from nearly 200 
typhoons within concentric, radial bands. They found that 

the most consistent averages relating storm movement and speed 
to surrounding winds lie within the 1-7 degree radial band. 

It appears that behond this zone, the influence of the shear 
in the environmental flow becomes more dominant. This 6-7 
degree "zone of influence" best delineates the cyclone/environ- 
mental flow characteristics. The component slope reversals 
Within this zone become better defined with height, with the 
Sharpest reversal at 250-mb. 

The component averages depicted in Fig. 6 were based on 
non-divergent wind components derived from the FNWC global 
band wind analysis. (Note that the hand-analyzed cases were 
Similarly treated.) The main difference between these fields 
and those from the hand-analyzed winds is the absence of the 
Slope reversal, except for the 50Q-mb u component. Only the 
850-mb components resemble the hand-analyzed cases. This 
may be attributed to FNWC's bogus initialization, which is 
most realisticat the lowest level, but weakens with height 
Meee eEisberry, 1977). Due to the higher reliability placed 
on the hand-analyzed cases (Fig. 5), along with similar re- 
sults from previous studies, all the area-averaged steering 
currents hereafter presented were calculated for a 7x 7 
(approximately 12 degrees longitude) square centered on the 


Srorm. 
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Peo UIMATLON OF THE EXPECTED MODEL STORM MOTION 
ieee initial: Observations 

Given the actual storm motion during the first 6 
hours of the forecast period, one must estimate the expected 
model storm track to determine the bias corrector that should 
be applied to the initial data fields. Running the model for 
approximately 6 hours would be ideal, but this would increase 
computer requirements. A plausible alternative may be to 
predict the model storm initial movement based on a number of 
prior case studies. Thus, a statistical regression equation 
was generated to estimate initial storm movement given the 
initial steering currents in the three layers. 

Development of the regression equation began with 
observing the initial 6-hour TCM vector of each test case. 

A 6-hour vector was chosen since maximum emphasis was placed 
upon improving the initial TCM storm movement. This vector 
had to be small enough to be correlated with the instantaneous 
steering current and yet large enough to represent the trend 
of initial storm movement. 

Hodographs depicting the average velocity vectors 
considered in the regression equation's formulation are shown 
in Figs. 7 and 8. One may compare the steering currents for 
each level, and the vertical mean current, with the TCM's 6- 
hour vector and the corresponding 6-hour best track vector. 
This latter vector is essentially the last 6-hour segment of 
the 12-hour bias input. All vectors were normalized with 


respect to the model storm vector which was prescribed to be 
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at 360°. Some of the results were Surprising. The first 
anomaly revealed by Figs. 7 and 8 is the lack of alignment 
between individual steering currents and the actual path of 
the storm. Comparison of the two figures shows that the 
500-mb currents have relatively similar relationships to the 
currents above and below. But the currents at this level 
were also most consistent with respect to individual cases. 
Comparison of the other levels revealed data differences. 

The reverse relationships of the 850- and 250-mb currents 
between the 1974 and 1975 cases were probably due to the use 
of hand-analyzed data in the 1974 cases, while the 1975 cases 
were derived from objectively-analyzed data. Specifically, 
the typhoon bogus applied in the 1975 cases resulted in rela- 
tively good flow fields at low levels, but poor fields aloft. 
The 1974 cases, on the other hand, exhibited superior upper 
level analyses due to the inclusion of aircraft and subjec- 
tively derived satellite data. 

The second characteristic common to Figs. 7 and 8 is 
that the model storm traveled about twice the speed of the 
averaged steering current. However, beyond the 6-h forecast 
periods, the TCM speed forecasts were too slow relative to 
the best track. The 1975 case steering currents were also 


ElleweCabout 2.6 eo 


),; which can be attributed to the objec- 
tive analyses. The vertical structures of the steering 
currents in the 1974 series appear to be more realistic. 


The individual steering current speed became greater with 


altitude and nearly matched the model's speed at 250-mb. 
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This observation suggests the importance of the upper level 
Steering currents for short-term motion forecasts and 
emphasizes the need for greater quantity and better quality 
Upper-level data. 

An alternative method of biasing was also investi- 
gated. .The recent study by George and Gray (1976) produced 
an empirical relation between the storm track and a steering 
current (mean flow) derived from rawinsondes. Their research 
included stratifications of the data according to storm lati- 
tude, translation speed and direction, intensity, and 
intensity changes with respect to the surrounding winds. 

The hodograph in Fig. 9 depicts their averaged empirical 
steering current relative to the storm's actual movement. 

It 1S important to note that the empirical steering vector 
was taken to be 86% of the mean 700-mb speed and lige toe che 
Paget OL the mean 500-mb direction. Notice the comparison 
between the averaged steering currents from Figs. 7 and 8 

and George and Gray's mean flow, which 1s essentially the 
storm steering current. This mean empirical vector and the 
1974 averaged current have nearly identical angular displace- 
ments with the latter being 30% slower. This close agreement 
was to be expected, since both vectors were constructed from 
raw data relationships, although the methods differed. The 
1975 objective analyses, however, produced an averaged current 
that had twice the angular displacement and was 33% slower 


than George and Gray's mean empirical vector. 
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STORM TRACK 





pawl.) EMpirrcal mean Flow relationship with 
respect to actual storm motion vector 
basecmonwGeoree fand oray (1975). Mean 
flow 1s 86% of 700-mb speed and i6° to 
the right of 500-mb direction, 
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2. Model Characteristics 

Experimental adjustments designed to alter the steer- 
ing currents by biasing the wind fields were first done with 
the hand-analyzed cases. Jit was found that the boundary con- 
ditions used for the TCM's initialization process (see Appen- 
dix), for the derivation of w and subsequent non-divergent 
wind fields, distinctly changed the resulting steering cur- 
rents. These adjustments altered significantly major portions 
e@emeene tlow fields. Comparison of Fig. 10 (original 850-mb 
wind field) and Fig. 11 (non-divergent 850-mb wind field) 
reveals the distinct boundary effects characteristic of this 
"channeled" model. New steering currents calculated after 
this initialization are presented in Figs. 12 and 13. Com- 
parison with the corresponding divergent steering currents 
(Figs. 10 and 1l, respectively) illustrates the changes. The 
effects of deriving the non-divergent components of the 1975 
cases were significant. Even larger vector changes occurred 
with the 1974 cases. The hand-analyzed fields contained 
Strong convergence/divergence areas, which were completely 
removed. An overall reduction of wind speeds was consistent 
with all cases studied. The dramatic differences between 
Figs. 8 and 13, and especially Figs. 7 and 12, dictate strong- 
ly that the steering currents must be calculated from non- 
divergent wind components to accurately reflect the storm/ 
environmental relationships within the model. 

Ss. Ihe Prediction Equations 
The purpose of the prediction equations as applied 


Pemeae wg 74 amd 19/5 eases was to eStimate the initial storm 
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Fig. 13. Same as Fig. 8 except with non-divergent winds. 





movement in the TCM from the initial steering current infor- 

mation. Statistical regression coefficients were calculated 

With a Bio-Medical Stepwise Regression program (Dixon, 1975). 
By contrast, the second adjustment method based on George and 
Gray (1976) came directly from the vector relationships pre- 

wemted in Fig. 9. 

UGMebeain EhewmacmumN 2nrormation from the individual 
steering currents, all the u and v components were used as 
independent variables (4a, b). The predicted components of 
the initial storm motion (Uy and Vu? were sums of their con- 
Stants (CU and CV, respectively) and the products of their 
Slx steering u, v current components and respective regression 
coefficients (A through F). Numbered subscripts denote model 


ayer. 


Uy, = CU + AUgcg + Begg + CUosg * DiVgsg + Eylsgq * FyYosg 89 CH? 
Vy = CV + AUgcg + Bue gg + OC Uo5g + OM esq t BYcgg t FyYo5g 86 SD? 


High confidence levels were found for the 1974 coefficients 
with the regression equations explaining 100% and 99% of the 
predicted u and v (respectively) storm components. The hand- 
analyzed data proved highly reliable even with the steering 
currents averaged from the non-divergent winds. The 250-mb 


v component alone explained 72% of the predicted V,, component. 


M 
It must be noted that only seven hand-analyzed cases were 
aeeelable | to produce these coefficients with such high confi- 


dence. More cases may have reduced their reliability. The 


Seneraenee "evels Of the 1975 coefficients were less than 
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those of the 1974 cases. These coefficients explained 78% 
and 68% of their respective Un and Vu variances. The de- 
creased effectiveness of the 1975 coefficients can partially 
be explained by the inferior wind fields initially supplied 
to the TCM, again emphasizing the need for an. improved 
cyclone bogus and the lack of better data for and analyses 
from the GBUA. In addition, the greater number of cases 
(41) explained by these coefficients would also account for 
their ineffectiveness relative to the seven 1974 cases. 

The two vectors depicted by Fig. 9 provided the neces- 
sary relationships needed to formulate an empirical equation 
based on the George and Gray results. This equation, how- 
ever, derives a mean flow from the bias input, not storm 
motion. The consolidated empirical vector, representing 
fractions of the 700-mb speed and 500-mb direction, was 
utilized as a vertically averaged steering current. This 
was necessary since the derived mean flow needed to be com- 
pared to the observed, vertically averaged steering current. 
A mean flow vector was then derived from the bias input, 
which represents the best track vector (Fig. 9), and the 
empirical relationships defining a corresponding storm mean 


imme (Sas b). 


g 


Empirical Mean Flow Speed = 86% x Bias Input 


Speed (5a) 
Empirical Mean Flow Direction = 16° + Bias Input 
Direction Goi» 


The bias input and its derived mean flow then reflected the 


identical relationships displayed by Fig. 9. 
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Pee ue BIAS CORRECTORS 
bts Pea Dproden 
A bias corrector is defined as the correction applied 
to a wind field to adjust the steering current to conform to 
the desired storm motion (Fig. 14a). For the 1974 and 1975 
storms, their initial motion was predicted from a regression 
equation. The bias correctors, Au and Av (defined Diya. niGhe “Oa. 


b), are 


Au (6a) 


Bias Input (u) - Uy 


Av = Bias Input (v) - V (6b) 


M 


the difference between their respective component bias inputs 
and predicted storm motion. These correctors were then 
applied to their respective component wind fields for steer- 
ing current adjustments as discussed in the next section. 
Ze Second Approach 

The relationships produced by George and Gray (Fig. 9) 
enabled derivation of a bias input's mean flow. Since a mean 
flow vector represents a vertical average, the derived mean 
flow and observed steering current components (U and V) were 
comparable. Thus, any differences between the two vectors 
(7a, b) represented the bias correctors needed to adjust the 


observed current to resemble the empirical current (Fig. l14b). 


Au = Empirical Mean Flow (u) - U (7a) 


Av = Empirical Mean Flow (v) - V (7b) 


An important difference between the derivation of the 


imj@esotCrOrebpias GCOormectors iS that the “storm track" derived 
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SC-STEERING CURRENT 
BI-BIAS INPUT 


BC-BIAS CORRECTOR 


PREDICTED ADJUSTED CURRENT 
STORM ___ 
TRACK 
(a) Sc 
EMPIRICAL CURRENT 
Bl 
(b) a 
meee lu. «6Sechematic of the bias corrector definition 


from (a) the first approach using the predicted 
initial storm motion and (b) second approach 
uSing the empirical steering current, 
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correctors (first approach) are differences between observed 
and predicted storm motion, whereas the "mean flow" derived 
correctors (second approach) are differences between the ob- 
served and empirical mean flow. Both sets of correctors, how- 
ever, adjusted the wind fields to produce a desired steering 


current average. 


Eepeelne BIASING SCHEME 
eee soltial Experiments 

Once the bias correctors were calculated, an effective 
method of applying them to the wind fields was investigated. 
Although the steering currents varied in the three layers, 
adequate information to differentially bias each level was 
not available. Therefore, the biasing scheme was simplified 
by distributing the bias correctors equally in each of the 
three layers. If a vertically dependent biasing scheme could 
be developed, improved forecasts may possibly result. 

First attempts at adding the derived bias correctors 
uniformly to the u, v fields over the entire area were not 
successful. Adding a uniform north-south component is incon- 
Sistent with the requirement for non-divergence in the initial 
Wind fields. It was therefore necessary to limit the biasing 
domain. The logical area would be that which originally 
defines the steering current, a 7 x 7 grid square. The im- 
provement in track forecasts obtained with the bias applied 
in a box-like function over the exact domain defining the 
Steering currents (Fig. 15a) was rather small. Several modi- 


fications were then tried by applying the maximum bias at the 
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Fig. 15. Examples of various horizontal distributions 
or the bias corrector with respect to the 
grid point nearest the storm center. In 
each case a unit value corresponds to the 
Jes ued bias COumeetor as am flg. 14, 
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storm center, and decreasing the magnitude to zero at three 
grid points from the center. Varied biasing functions from 
linear to cosine (Figs. 15b, c, and d) were tried in an effort 
to eliminate the discontinuities believed to have adversely 
affected the wind fields by the box-like function. [In addi- 
tion, the magnitudes of the corrections were systematically 
varied with each function tried. These efforts proved un- 
fruitful as many of the forecast tracks were ill-behaved, and 
in some cases, the tracking routine could not follow the 
center. Analysis of the forecast wind fields showed that the 
bias had altered the vortex so drastically that its basic 
structure was not maintained during the initialization 
processes. 

The two-fold objective of adjusting the steering cur- 
rent and retaining the vortex indicated that the modification 
should be concentrated in the "zone of influence". The zone 
of influence was previously defined as the region about a 
vortex which is affected by both the storm and the environ- 
ment. Application of the maximum bias just beyond the zone 
of influence (Fig. 1l5e) proved most effective. Many varia- 
tions in placement of the maximum See finally resulted in the 
form shown in Fig. 15e, where the maximum bias was at the 
fourth grid square. Although the correction values applied 
to this function were later modified, the placement of maximum 
bias was kept fixed. As an illustration of the effects of the 
modified biasing form used (described in Results section), 
the response of the wind fields to this type of biasing is 


apparent upon comparison of the unbiased 850-mb winds in 
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Mego i) and the biased winds in Fig. 16 for the 4 August 1974 
case. The major effect was the alteration of the system's 
zone of influence. The winds near the vortex center under- 
went minor change, but the flow region (zone of influence) 
imeeracting with the surrounding flow resulted in maximum 
change. In this case, the surrounding flow patterns were 
altered significantly. The biasing literately reformed the 
adjacent pressure systems to insure the storm's northward 
movement. This was the result of the substantial bias applied, 
combined with the restrictions imposed by the model's boundary 
conditions. Although seemingly excessive, these biased wind 
Fields resulted in significant improvement in the track fore- 
cast, as will be discussed in detail in the Results section. 
Because the bias introduced both vorticity and diver- 
gence in the initial wind fields, a second streamfunction 
solution must be calculated to generate non-divergent wind 
components. However, this result deviated from the desired 
Steering currents. Comparison of the biased and subsequent 
non-divergent fields revealed wind speed reductions as great 
as 50% in regions of maximum bias. This phase of the TCM 
initialization consistently decreased the magnitude of the 
bias by blending sharp gradients and eliminating convergent/ 
divergent patterns introduced by the bias. Therefore, to 
counteract these effects, the biasing scheme was applied 
twice. That is, after the original initialization and calcu- 
iielonm Of the bias correctors, biasing and reinitialization 


was repeated. This method improved the forecasts, whereas 
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additional applications of the bias altered significantly 
the winds and did not result in consistent improvements. 
2. Differential Adjustments 

The procedures performed in this section were purely 
developmental. The basic method was trial and error. These 
steps were essential since the characteristics of each TCM 
(1974 versus 1975) are highly dependent on the latitudinal 
domain. This is a direct consequence of the no-flux condi- 
tions imposed on the north-south boundaries. The TCM domain 
is not relocatable during the prognostic phase. In some 
cases, the channel walls restricted northward movement as 
Seposed tO westward movement. The average model storm track 
in the 1974 cases without biasing was 269°, while the average 
best track direction was 296°. Corresponding values for the 
1975 cases were 293° and 338°. Thus, one of the effects of 
the bias adjustment is to offset this model-related hindrance 
to the north-south storm motion. 

The method employed to overcome these north-south 
restraints was differential component corrections. Beginning 
Poem the basic biasing function of Fig. 15e, the amplitude 
was systematically varied. Each 48-hour test forecast was 
plotted against the unbiased track and the best track. The 
ultimate goal of the biasing was to obtain an initial six-hour 
track nearly parallel to the best track without producing an 
unrealistic 48-hour track. At the minimum, the biasing 
should result in a track that matched the unbiased forecast. 

The first bias application, as in Fig. 15e, used the 


original bias correctors Au and Av for the maximum values, 
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or the "peak" values. Then, holding the peak u-component con- 
stant, the peak v-component was either increased or decreased 
until the best solution was observed for the seven cases 
tested. Finally, the new peak v-component was held fixed as 
the peak u-component was varied. This process was repeated 
Several times until no further improvements were observed. 
During the development of the biasing function illus- 
trated in Fig. 15e, it became apparent that the v-component 
bias would have to be much larger than the u-component to 
offset the TCM's north-south boundary conditions. The rela- 
tive magnitudes and peak values of the u, v bias correctors 
were different for the 1974 and 1975 cases, for each used a 
different domain size. It might be anticipated that the 
replacement of the channel boundary conditions with a relocat- 
able model with open boundaries would eliminate some of these 
problems. There should also be little or no need for 


differential component biasing. 
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Ve REsULIS 


PweeeiHe 1974 CASES 

The main purpose in tests with the 1974 cases was for 
preliminary development of the biasing scheme. Results from 
these cases will be presented on a relative improvement 
basis only. Although there were small errors associated 
With the initial positions, they were subtracted from all 
results presented. However, these errors will be included 
With the 1975 results as operational feasibility is examined. 

The 1974 cases provided an excellent sample because the 
data were hand-analyzed and the statistical regression co- 
efficients resulted in very accurate estimates of the initial 
model storm motion. These cases were therefore used for pre- 
liminary testing which finally resulted in the basic biasing 
Function of Fig. 1l5e, as described above. An example of 
the results thus far is shown in Fig. 17. Because the 
original unbiased forecast was so poor, this case (Typhoon 
Gilda, 04 August) was of special interest. These TCM fore- 
casts were plotted every 6 hours to illustrate the temporal 
stability within the model and the short-term effects of the 
various biasing functions. The unbiased 48-hour track is 
labeled as B in Fig. 17, and may be compared to the best 
track (A). The first attempts with the box-like, cosine and 
basic biasing functions are shown in tracks C, D and E. The 


maximum values applied to the functions were the unit Au, Av 
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Fig. 17. Model forecast tracks from 00 GMT O04 Aug 
1974 versus (A) best track and (B) unbiased 
and for a (C) Box=-like, (D) Cosine and 
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bias correctors. The box-like biasing (C) resulted in an 
abnormal track forecast with an average speed of about one 
aS pe Although the cosine function (D) generally improved 
the forecast, increased weighting of the bias distorted 
seriously the cyclone wind fields. The forecast resulting 
from application of the basic biasing function (E) was not 
readily apparent as the optimum choice. However, upon in- 
creasing the weighting of the bias correctors, this was the 
only function which remained stable, even when applied with 
peak values an order of magnitude greater than the unit bias 
correctors. These corrections were easily tolerated by the 
TCM, since the solution for non-divergent wind components 
during reinitialization quickly reduced the added bias. 
That 18, the wind fields needed to be over-adjusted to make 
the corrections effective. 

The seven 1974 cases were used for testing the basic 
Pacing function through differential component modifications. 
Analysis of track E (Fig. 17) suggests the need for a stronger 
northerly bias. The unit biaS correctors for this case were 
Au = 0.8 ms + and Av = 3.8 ms-*. Not as obvious was the need 
for less easterly bias since Au seemed to have reduced the 
speed of track E. Therefore, to force greater northward move- 
ment while maintaining acceptable forecasts from all test 
cases, the Av peak biasS was increased as the Au peak bias 
was decreased. The resulting biasing function for the v com- 
Mement 1S shown in Fig. 18. It was found for this limited 


Sample that no u bias or negative v bias was necessary. The 
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general shape of this curve resembles the basic biasing func- 
tion of Fig. 15e, but the slightly non-linear slope was found 
to produce better results. Because of the large v bias used, 
this was a successful effort to reduce distortion of the 
central vortex by the bias. As shown by Fig. 18, the v-com- 
ponent corrections for a Av of 4 mst would be extremely large. 
Even a bias corrector of 2 ms~2 would require a peak value of 
yo ms7*. It should be emphasized that the very large forcing 
of the v-component was peculiar to the small domain version 
e-eene TCM used. 

The final results for all four Gilda cases are shown in 
Fig. 19. The actual storm path 1s plotted for the 0000 GMT 
positions. Note that the storm actually recurved while all 
of the unbiased forecast tracks were basically westward. Each 
biased forecast track is a result of biasing twice using 
weighting function represented in Fig. 15e-. These biased 
forecasts resulted in little and moderate improvements for 
days 01 and 02 respectively. Excellent results were obtained 
for both days 03 and 04. The significant accomplishment of 
this biasing scheme was to make the storm recurve. The re- 
sults for the other three cases (not plotted), for Typhoon 
[emewon 25, 26, and 27 November 1974, did not significantly 
alter the unbiased forecasts. The biased forecasts for Irma 
had an average error increase of only 37 km for the 48-hour 
period. Since the actual storm and unbiased tracks both 


moved westward, it was difficult to improve forecasts that 


were already very good (Ley, 1975). 
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Fig, 19. Model forecasts with bias distribution as 


in Fig. i5e for the 01-04 Aug. 1974 cases, 


62 





Table I presents the forecast error averages for both this 
biasing method and the empirical approach based on George and 
Gray (1976) as described above. The number of biased fore- 
casts which equalled or exceeded the unbiased forecasts are 
also tabulated for each six-hourly forecast interval. Figs. 
20 and 21 depict graphically the averaged errors and standard 
deviations (respectively) from Table I. As originally required 
of this biasing scheme, the initial forecasts did improve. [In 
fact, these cases resulted in track forecast error decreases 
of about 40% during the first 30 hours and about 25% at 48 
hours. The method based on George and Gray (1976) decreased 
forecast errors by about 20% on the average. The official 
JTWC 24- and 48-hour forecast errors for the same seven fore- 
cast periods are also plotted in Fig. 20 for comparison. Not 
only did biasing by the first approach significantly decrease 
the TCM forecast errors, but also gained more accuracy over 
JTWC's official forecasts. Of equal significance, if not more 
so, are the standard deviation curves in Fig. 21. This graph 
shows that by biasing, especially by the first approach, 
there was less than half the variance of that resulting from 
the unbiased forecasts. The biasing resulted in tracks that 


were much more stable during the entire forecast period. 


pees 1975 CASES 

These 41 (originally) cases were initialized from the 
FNWC objectively-analyzed data (GBUA). All cases were used 
for calculation of the prediction equation regression coeffi- 


Clients. Of the original series, three cases were not 
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TABLE IL 


Method Unbiased Pits eras SeGend Jbras 

Spearistic at al ied I na 

Cee hr C/) 43 Z2 (5) 26 (4) 
meant C7) 79 44 Ss) 62 (4) 
mom 4 /) 119 oD oe 85 CS) 
Zope (7) 162 oa CS.) ae (4) 
coenr C7) DANS Ze is) 164 Ca) 
Seemr C/) ZO 189 (2) Z0F Cy) 
Mawr (/) 5 255 (5) Zi Co) 
48 hr (7) 426 724 co) 394 (4) 


Mean 6-hourly forecast errors (km) (I) based 

On Best track for the 1974 cases for the 
unbiased, first bias and second bias approaches. 
The error at each 6-h interval is normalized by 
the initial position error. Also tabulated 

are the number of cases (II) which equal or 
exceed the unbiased error. 
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forecastable because the TCM's time step was not adjusted 
for storms occurring in the higher latitudes. Notable differ- 
ences between the two series of tests (1974 and 1975) included 
the regression coefficients, the final forms of the biasing 
Funetions, the use of operational data versus hand-analyzed 
data, and model domain sizes. The results from an entire 
season of tropical storms should reveal how flexible this 
biasing scheme is. 

The same basic biasing function from the 1974 cases (Fig. 
lose) was used for obtaining the final biasing form best 
Suited for the TCM version at FNWC. Again, differential bias- 
ing experiments were conducted by systematically adjusting 
the peak value of one bias corrector while holding the other 
fixed, and then vice versa. Because the FNWC TCM was four 
grid points wider and longer, a u-component bias was neces- 
sary, as had been expected. But after application of the 
resulting biasing function, the error showed less improvement 
than expected. In fact, the blasing increased the average 
forecast error, but only about 5%. Some of the error may be 
attributed to the estimate of the initial model storm move- 
ment, Since the statistical regression equation based on the 
1975 data only explained about 70% of the variances, as com- 
pared to 99% for the 1974 cases. However, the main source of 
error appears to be the initial data fields used by the TCM. 
First, several of the storms were not properly bogused. This 
was discovered by comparing the JTWC operational surface 


analyses with the initial 850-mb wind fields. These same 
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wind fields also revealed discrepancies within the original 
GBUA analyses (see Elsberry, 1977). Another possible cause 
was that several cases with maximum winds below 50 kts were 
included in the sample. The operational version of the TCM 
is only used for storms exceeding 50 kts. One might also 
note that the biasing scheme had little beneficial effect for 
the 01 August 1974 case in Fig. 19, when Gilda had 40 kt 
maximum winds. Analysis of the individual 1975 cases sug- 
gested that the biasing function would be more effective if 
the peak biasing value was proportional to storm intensity. 
However, in the form tested here, the biasing did not depend 
on storm intensity. During these tests, it was observed that 
a Significant division of results from weak versus strong 
biasing occurred about the 60 kt maximum wind speed. There- 
fore, in an effort to test this biasing scheme on a more-or- 
less congruous sample, the only cases (24) tested were storms 
having maximum winds of 60 kts or greater and those with good 
initial analyses and vortex boguses. Since greater importance 
was Dplaced on the more developed storms, the biasing function 
was adjusted to perform best with the stronger storms. 

Having thus defined the test sample, differential bias 
testing produced the biasing function represented in Fig. 22. 
Although these curves show the bias that would be applied if 


both bias correctors, Au and Av, were 4 mst, most of the bias 


Serrectors were about 1-5 ms7*. The negative correction 
values are also shown. Note that the slopes are identical 


to that of the basic biasing function. This sample (linear) 
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function was sufficient, since the large bias required for 
the 1974 cases Cand resulting non-linear slope) was unneces- 
sary due to the wider TCM channel. In general, only half as 
much v-component bias was required in this TCM as compared to 
that needed by the 1974 series TCM. The positive (eastward) 
and negative (westward) u-component corrections required for 


1 ; , 
uUn1lt corrector are only one-third and one- 


the same 4 ms 
twelfth the strength needed by an equivalent v-component unit 
Serrector. These functions indicate that a combination of 
mostly northerly and some easterly bias will aid storm recur- 
vature in the model. Fig. 23 illustrates graphically the 
respective weights of these differential corrections. The 
central x-y axes represent the unit bias corrector values, to 
be applied at the fourth grid square as in Fig. 22. By 
imagining Fig. 23 relative to a channeled TCM, one can visu- 
alize the greater bias applied by v-component and positive 
u-component correctors. Note that the peak u biasing scales 
are four times smaller than the v biasing scale. The negative 
v-component was biased just as heavily as the positive compo- 
nent. It was also observed that large biasing by one compo- 
nent relative to the other component often resulted in poorer 
forecasts. This was especially true if this other component 
was less than l ome Therefore, the dashed lines determine 

a component's peak value under conditions determined by the 
Other component's unit value. In general, if Au was less 

than % eae then the peak v value would be reduced by 75%. 


If the absolute value of Av was less than 1 ms7-, then the 


peak u value would also be reduced by 75%. 
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The error statistics resulting from the 1975 test cases 
are presented in Table II. The unbiased and biased TCM and 
the JTWC 48-hour forecast errors are shown in Fig. 24. The 
biasing scheme attained consistent improvements over the un- 
biased cases throughout the forecast period. Forecast errors 
with the biased model versus the unbiased model were decreased 
about 30% during the first 30 hours and about 15% after 48 
hours. The JTWC official forecasts were equalled or exceeded 
by the TCM forecasts with biasing. Average standard devia- 
tions for the unbiased and biased forecasts are shown in 
Fig. 25. The biasing decreased storm track variations by 
about 45% for the first 36 hours. As with the 1974 case re- 
sults, the decrease of the standard deviation error was the 
most significant achievement of biasing. 

The above comparison of biased and unbiased forecasts 
was in terms of the relative errors. That is, the positional 
errors originating during initialization were not included. 

A comparison of the same forecast schemes as in Fig. 24, 

with the inclusion of the errors due to the poor initial 
werus Of the storm location, 1S shown in Fig. 26. Even with 
this crude model, the bias adjustments produced results after 
30 hours nearly equivalent to JTWC official forecasts. In- 
cluding the 41 unbiased cases, the average initial position 
error for the TCM was 85 km, whereas the official JTWC 
position error was only 34 km. Therefore, the effects of 
biasing cannot be expected to improve this error, since the 


blasing procedures are based on the wind fields provided by 
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Method Unbiased FifSsieolas 

Statistic T i 0s 
Q6 hr (23) 8 | 55 cae 
12 hr (23) 128 69 (21) 
18 hr (23) 182 125 a 
24 hr (22) DD 17 is 
30 hr (18) 289 214 (13) 
36 hr (17) 345 267 (13) 
42 fr (13) 349 311 ( 8) 
48 hr (13) 413 363 (10) 


came aS Table I except for the 1975 cases and 


without the second bias method, 
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the initial storm bogusing. But it should be noted that the 
biasing scheme did introduce an additional initial error of 
about 15 km. As a result of determining the average steering 
current and applying the bias corrections based on the grid 
point nearest the storm center, the adjusted steering current 
may have become asymmetrical. Adding too much bias in one 
quadrant evidently alters the storm circulation center. This 
problem should be solved by using interpolation schemes for 
both the steering current averages and the application of 
bias corrections. 

Of the 24 test cases, two of the forecasts with the bias 
did not complete the 48-hour period. One case could not be 
tracked after 24 hours. No apparent explanation was found, 
but the bias may have adversely affected the vortex structure. 
The second case was comparatively weak (70 kts), which suggests 
that an intensity~proportional biasing function may have 
reduced the bias applied and possibly permitted a 48-hour 


merecast. 
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VI. CONCLUSIONS 


A primitive-equation tropical cyclone forecast model (Ley 
and Elsberry, 1976) is being tested at FNWC (Hinsman, 1977) 
With operational data. The objective of this research was to 
improve the initial model storm tracks by incorporating the 
observations of actual storm motion that would be available 
before the model is run. Development consisted of testing 
two series of initialization data. The 1974 data cases were 
hand-analyzed by JTWC, while the 1975 data cases were objec- 
tively analyzed by FNWC. Two biasing approaches were 
developed. The first was based on a statistical regression 
prediction of initial storm motion while the second used an 
empirical relationship (see George and Gray, 1976) to define 
the steering current. 

The largest improvements occurred within the 12- to 36-h 
periods, and the effects continued beyond 48 h. This result 
was found with both the 1974 and 1975 cases. The relatively 
large errors for the first 6-h forecasts during 1975 can be 
eeeeer buted to two effects. First, most dynamical models 
Pomme an adjustment period of 6 to 12 hours. Second, it 
requires some time for the effect of the bias adjustments at 
relatively large distances from the storm center to take 
effect. Some attempts were made to apply enough bias to make 
the initial 6-h forecast nearly parallel to the desired 


track. However, such extreme biasing always resulted in an 
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abnormal track during later periods because excessive bias- 
ing of the initial wind fields resulted in severe deforma- 
tion of the synoptic features. 

It should be pointed out that the weakest link in this 
biasing scheme is the prediction of the initial storm motion. 
The regression coefficients have not been tested with an in- 
dependent sample of storms. Application of the statistical 
equations in the future may not produce comparable results. 
Even the estimation of the steering current based on George 
and Gray (1976) is based on historical correlations and may 
not be relevant during an abnormal storm season, such as 
during 1975. A superior approach would be to run the model 
for six hours, which can then be compared to the actual storm 
motion. This procedure would require some additional com- 
puter resources but would Premina te the uncertainty asso- 
Clated with the statistical regression equation. These 
regression coefficients may also be geographically dependent 
and thus may be limited to the western Pacific region. 
Elimination of the prediction equation would therefore permit 
application of this biasing scheme to other regions. 

A major source of error ‘encountered with the 1975 test 
cases was the quality of the data fields originally supplied 
to the model (see Elsberry, 1977). The objectively-analyzed 
wind fields (GBUA) were significantly inferior to the hand- 
analyzed cases. The storm bogus, which initially positions 
and defines the vortex, not only misplaced the storm center, 


but the intensity and vertical structure were unrealistic at 
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times. Further attempts to improve this biasing scheme 
Should not be attempted until the bogusing routine is improved. 

This thesis found an effective method for applying the 
bias horizontally, but there is still some uncertainty as to 
how to apply the bias in the vertical. It is expected that 
a vertically dependent biasing scheme will produce better re- 
Sults. This should be especially applicable to the 250-mb 
level where cyclonic outflow patterns may need special treat- 
ment. It is also suggested that a simple intensity-propor- 
tional biasing function should produce better results, because 
storms with less than 60 kts were more sensitive to strong 
bias corrections. 

This biasing scheme can improve the dynamical model fore- 
cast accuracy by using the actual storm motion, which is 
currently available when the operational model is begun. The 
potential abilities of this method were observed with the 1974 
series results. Operational use is presently possible if the 
most recent storm track 1s well known. Elimination of the 
prediction equation and the inclusion of a storm intensity- 
weighting biasing function should lead to improved forecast 


guidance from tne dynamical model. 
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PaO ene Or CAL CYCLONE MODEL 


fee LHe MODEL 

The primitive equation model developed by Elsberry and 
Harrison (1971) and Harrison (1973) was used as the basis 
for these experiments. Although this model is capable of 
triply-nested operation, insufficient time limited testing 
to the uniform, coarse-mesh grid version. The model equa- 


mlons are: 











2) 0 KO) OT see OT 5 
ser cL(u) + fv - Mh + pe (A-1) 
ee = -L(v) - fu - M oS Txy + ud (A-2) 
a = -L(@) (A-3) 
pees 5 ESOT O02 aaa 
= = -M*(S. (2) + = (41 (as) 
se = @ C, 5 (Pony RIC (A-6) 
where 
L(S) = Mote. (44) + = (YS) + x5 (ws) 


L(S) represents the flux divergence of any scalar quantity 
S. The other meteorological symbols used above can be found 


imeche “List of Symbols". 


oe 





The linear computational stability criterion for two- 


dimensional equations governing simple wave motion 1s 


Carr 


ee 5 0e (A-7) 
where 
eee oe piece clOci ty Cf the fastest eravity wave 
At = time increment 
Ax = horizontal grid increment 


Eq. (A-7) dictates a maximum time step of 480 seconds for 

this model. An increase to 800 seconds was achieved by time 
averaging the pressure gradient term of the momentum equations 
(see Ley, 1975). The 800-second time step was used for test- 
ing the 1974 cases, while a 600-second time step was used for 
testing the 1975 cases as the relocatable grid was extended 
northward where Ax 1s reduced. 

The initial step was forward in time. All subsequent 
iterations used the leap-frog scheme. Parameterization of 
latent heat release was not possible because moisture fields 
were not available. Consequently, friction was also neglected, 
so that the TCM forecasts storm movement primarily by 


advective processes. 


fee lit GRID 
The uniform, coarse-mesh grid 1S a mercator projection 
ene at 22-5 N. The horizintal grid increment was nominally 


220 km. The grid extended 28 grid points east-west and 20 
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grid points north-south for the 1974 cases. An extra four 
grid points east-west and north-south were used for the 

1975 cases. The grid was oriented such that each storm was 
initially centered near the 17th column and the 10th row to 
allow for movement into the NW quadrant. Fig. A-1 shows 

the vertical distribution of dependent parameters. Although 
the variables are staggered in the vertical, horizontal 


Space staggering was not utilized in these experiments. 


meee SOUNDARY CONDITIONS 
Boundary conditions on the northern and southern walls 
were the no-flux type after Elsberry and Harrison (1971). 


imrese conditions are 


real = >. (A-~8a) 
M. 
SO ae’ ie 
aa} = M we (A-8b) 
Use] = Us (A-8c) 


Such a representation actually places the wall between the 
outer two grid rows on the poleward and equatorial sides. 
Mass 1S conserved along these boundaries since the v compo- 
nent is always zero at the wall; in essence, the wall is a 
streamline. Thus, the 18 grid rows of extracted data were 
expanded to 20 rows with the inclusion of the no-flux 
boundaries. 

The east-west boundaries were made cyclic after 


Krishnamurti (1969). This procedure, applicable to the 
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Fig. A-1. Vertical distribution of dependent variables 
and pressure levels for the three-dimensional 
medel -Caftter Harrison, 1973)4 
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1000-mb height field and the 850-, 500-, and 250-mb u and v 
fields, adds a buffer zone in the east-west direction to 
absorb the impact of forced continuity. For the 1974 cases, 
the original data of 23 longitudinal columns was expanded to 
28 columns by adding one column to the west side and four 
columns to the east side. The values in this buffer zone 


were determined by 


08, = > 944 (A-9a) 


Soe SO 


Cac 
ied la 


and then linearly interpolated to fill values in columns 25 
through 27. The grid used for the 1975 cases was expanded in 
a Similar manner. Here S represents any of the original 

data fields mentioned earlier. 

The vertical velocity at the upper boundary is equal to 
Zere and is calculated at levels 5, 3, and 1, Fig. A-l, 
through downward integration of Eq. (A-5). Optimum results 
were obtained when, in Eq. (A-4), w = was replaced by the 
equivalent expression -wa. Further solution improvement was 


attained when w in Eq. (A-4) was represented by an averaged 


Veaene Trom 


ee ea eee ae 3278 SAH TOD 


DeeeeENWC'S TROPECAL CYCLONE MODEL 
The 1975 storms were run on a semi-operational basis us- 


ing FNWC's TCM. This model was modified by R. Perry of the 
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Naval Environmental Prediction Research Facility (NEPRF) and 


Lt. D. Hinsman, USN, of FNWC to run with operational data. 


It has two significant improvements over the version used to 


forecast the 1974 storms. 
filtering and a direct solver 
relaxation (SOR) (see Rosmond 
the Poisson equations for the 
model is currently supporting 


(Hinsman, 1977). 


Meee DIAGNOSTIC PHASE 


The model incorporates Robert time 


instead of successive over- 
and Faulkner, 1976) to solve 
ie eZ a On! Pp rOC@es > ae l las 


JTWC on an operational basis 


The relative vorticity fields were obtained from the ob- 


served u and v components 


ie BU Ss 


V Pe) u 
mM - by SH 


)] (A-11) 


Subsequently the streamfunction, wt, was found by sequential 


over-relaxation of the expression 


Belmtion of Eq. 


values on the northern and southern peripheries. 


(A-12) 


(A-13) 


(A-12) requires specification of boundary 


Lo remann 


consistent with the no-flux boundaries, described by Eqs. 


(A-8b and A-8c), constant streamfunction values were defined 
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for the northernmost and southernmost two rows. 


Elsberry and Harrison (1971), 


Following 


zero was chosen for the 


Southern rows, and the northern values were calculated from 


y = 


Orr hh 


Here u 


south © 


Gl 


(=P) Ay (A-14) 


W 1s the mean zonal wind component averaged over the 


entire grid and Ay is the distance between the north and south 


walls. 


An appropriate balance equation can be derived through 


partial differentiation of Eq. 


Eq. (A-2) with respect to y. 


(A-1) with respect to x and 


Addition of these two equations, 


plus rearrangement of terms, leads to 








lui LEG) V U 
Dae v Ja 5 v ee a dW 
Z 2 
Z. 0 ie Zen CKO) 
—Mu, — () —M s+ SS) (A-15) 


To begin the prognostic stage 
@mmeaq. (A-15) is set equal to 


were neglected. Manipulation 


Poisson equation, 


Weak 
oy M 


2 5 Cae) 


7 v 
ae a =a) 


smoothly, the left hand side 


ZOrO tame Sir aeeradistaiisnive vet fects 
of Eq. (A-15) leads toa 
“by Yb eae 
+ fio. ue soo wy, 3557 
(A-16) 


that is readily solved with sequential over-relaxation or 


Other direct methods. 


In this case, the mean geopotential 


value on the south wall was derived hydrostatically (Eq. 
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(A-6)) based upon a mean climatological temperature sounding. 
An average geopotential value for the north wall was found 


through integration of the v momentum equation (noting that 


OV 


coin 0) over the region. 





Se Te Oe Yo Ay (A-17) 


Again, Ay is thenorth-south grid distance. To assure mainten- 
ance of the no-flux boundary conditions, Eq. (A-8a) was now 
eepliecd at all levels. Having found the balanced geopotential 
fields from the upper level wind fields, only the 1000-mb 
geopotential field and the temperature fields at levels l, 2, 
4, and 6, Fig. A-1l, need to be specified to begin the prognos- 
tic phase. Since the initial wind fields are non-divergent, 
the vertical velocity 1s everywhere zero. 

The 1000-mb geopotential field was obtained hydrostatical- 
ly from the 850-mb values. This requires that the 1000-mb 
and 850-mb potential temperature fields are defined. From Eq. 
(A-6), recall that the thickness between two levels in the 
atmosphere is proportional to the mean potential temperature 
of the layer. Constant potential temperatures for the 1000- 
mb north and south walls were obtained from climatology. A 
linear gradient at 1000 mb was forced to fit the boundary 
values. It was assumed that the 850-mb potential temperatures 
were everywhere 4°K warmer than the corresponding 1000-mb 
values. This assumption resulted in a linear meridional gra- 


dient in the 1000/850-mb layer mean potential temperature 
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field. Finally, the 1000-mb geopotential field and the 
200-mb and 250-mb potential temperature distributions were 


obtained using Eq. (A-6). 
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